! ***************************************
finish
/clear
/filn, e344
/title, Stress analysis of skirt supporting zone of hydrogenation reactor
/units,si
/dev,font,1,Courier*New,400,0,-13,0,0,,,
/dev,font,2,Courier*New,400,0,-14,0,0,,,
! ********* 趨*********                                        
pi=3.1415926                            ! Բ
tongt_di=2813                           ! Ͳھ
tongt_t=87                              ! Ͳں
tongq_h=615                             ! Ͳ嵽ȹ֧дĸ߶
fengt_t=52                              ! ͷں
fengt_di=2833                           ! ͷھ 
xyhr1=20                                ! ڲԲ뾶
xyhr2=20                                ! Բ뾶
h1=1500                                 ! ȹ߶
h2=1000                                 ! Ͳ峤
xqunz_h=70                              ! hζͼȹӲֱ߸߶
qunz_t=22                               ! ȹں
yqunz_h=690                             ! ͷȹ䱣²λ                                     
fengt_do=fengt_di+2*fengt_t             ! ͷ⾶
baowen_t=180                            ! ²
Ktt=35*1e-3                             ! ͲĴϵ
Kbw=0.134*1e-3                          ! ²㴫ϵ
Hair=12*1e-6                            ! ϵ
Hyt=1000*1e-6                           ! ϶ϵ
Tair=20                                 ! ¶
Tyt=347                                 ! ¶
                                       
!****************ǰ***************************
/prep7
ET,1,PLANE77                            ! 嵥ԪΪ8ıȵԪ
KEYOPT,1,3,1                            ! ģΪԳƽṹ
MP,KXX,1,Ktt                            ! ͲĴϵ
MP,KXX,2,Kbw                            ! 屣²Ĵϵ
!****************ģ***************************
k,1,0,0                                 ! ĵ
k,2,tongt_di/2+tongt_t,0                ! hζͼ϶
k,3,tongt_di/2,0                        ! hζͼ϶ڲ
k,5,tongt_di/2+tongt_t,-tongq_h         ! hζͼȹӲ
k,6,tongt_di/2+tongt_t-qunz_t,-tongq_h	! hζͼȹӲڲ
k,7,kx(6)-xyhr1,-(tongq_h-xqunz_h)	! hζͼȹӲԲ

l,2,3                                   ! 
l,2,5
l,5,6
circle,7,xyhr1                          ! ȹӲԲԲ
ldele,5,7,1,1                           ! ɾԲ

l,6,4                                   ! 
k,7,0,ky(8)                             ! ¶ؼ7ΪԲȸߵĵ
l,7,8                                   ! hζͼԲߵĸ
circle,1,fengt_do/2                     ! ڶӦԲ
ldele,7,9,1,1                           ! ɾԲ

lsbl,10,6                               ! ߶
ldele,7,,,1                             ! ɾ߶

l,8,10                                  ! ٴԲߵĸ
k,30,kx(10)-xyhr2,ky(8)-xyhr2           ! 帨ؼ
k,40,kx(8),ky(30)                       ! 帨ؼ
 	
l,40,30                                 ! ɸȷڲԲ
circle,1,fengt_do/2+xyhr2               ! ԲȷڲԲ
lsbl,7,12                               ! 󽻵ȷڲԲ
ldele,13,14                             ! ɾ߶
kdele,30                                ! ɾĹؼ
ldele,9,11,1,1                          ! ɾ߶

circle,14,xyhr2                         ! ڲԲԲ
ldele,7,,,1                             ! ɾĹؼ߶
ldele,10,11,,1                          ! ɾĹؼ߶

l,1,14                                  ! ɸ
l,9,14                                  ! ɸ
k,5000,kx(14),ky(14)                    ! ٶڲԲ
lsbl,9,7                                ! ߶
ldele,12,,,1                            ! ɾ߶

l,7,14                                  ! ɸ
lsbl,6,10                               ! ߶
ldele,12,,,1                            ! ɾ߶
lsbl,8,7                                ! ߶
ldele,6,,,1                             ! ɾ߶

k,1,0,0                                 ! ¶
circle,1,fengt_di/2                     ! ڱڵԲ
ldele,6,8,1,1                           ! ɾ߶

k,40,kx(3),ky(7)                        ! 帨ؼ
l,3,40                                  ! ɸ
l,12,14                                 ! ɸ
lsbl,12,6                               ! ߶
ldele,8,,,1                             ! ɾ߶

ldele,3                                 ! ɾ߶
k,100,kx(5),ky(5)-h1                    ! ȹ߶εĹؼ
k,200,kx(6),ky(6)-h1                    ! ȹ߶εĹؼ
l,100,200                               ! ȹ߶
k,300,kx(3),ky(3)+h2                    ! Ͳ߶εĹؼ
k,400,kx(2),ky(2)+h2                    ! Ͳ߶εĹؼ
ldele,1                                 ! ɾ߶
l,300,400                               ! Ͳ߶ 
l,300,11                                ! Ͳڱ߶
l,400,2                                 ! Ͳ߶
l,5,100                                 ! ȹ߶
l,6,200                                 ! ȹڱ߶
k,800,kx(2),ky(11)                      ! Ӧ·ؼ
k,900,kx(5),ky(4)                       ! Ӧ·ؼ
k,6000,kx(300),ky(300)-400              ! Ӧ·ؼ
k,7000,kx(400),ky(6000)                 ! Ӧ·ؼ
l,1,5000                                ! ɸ
lsbl,13,15                              ! ߶
ldele,2                                 ! ɾ߶   
a,12,14,10,11,300,400,2,800,900,5,100,200,6,4,8,9,7     ! 

k,8000,kx(400)+baowen_t,ky(400)         ! Ͳλñ²ؼ
l,400,8000                              ! Ͳλñ²߶
lsel,s,loc,x,kx(400)                    ! ѡͲȹ߶
adrag,all,,,,,,18                       ! Ͳȹı²

k,8100,kx(8000)+100,ky(8)-yqunz_h       ! ͷȹ䱣²λø
k,8200,0,ky(8)-yqunz_h                  ! ͷȹ䱣²λø
l,8100,8200                             ! ͷȹ䱣²λø
 
k,8300,0,0                              ! ¶
circle,8300,fengt_do/2+baowen_t         ! ɰͷ²߶εԲ
ldele,31,33,1,1                         ! ɾ߶

lsel,s,,,12,14,2                        ! ѡ񱻷ָ߶
lsel,a,,,28,34,6                        ! ѡ񱻷ָ߶
lsel,a,,,30                             ! ѡָ߶
lsbl,all,30                             ! ָ߶
ldele,36,,,1                            ! ɾ߶
l,24,21                                 ! ɷͷȹ䱣²λ߶
alls
l,22,12                                 !  
a,12,7,9,8,4,6,21,24,22                 ! ɷͷȹ䱣² 

aglue,all                               ! ճ
aesize,all,10                           ! 嵥Ԫߴ
mshkey,0                                ! ʷ
asel,s,,,1                              ! ѡṹ
mat,1                                   ! 趨
amesh,1                                 ! ʷ
lrefine,4,,5                            ! Բ
lrefine,11,,5                           ! ڲԲ

alls
asel,u,,,1                              ! ȥṹ
mat,2                                   ! 趨
amesh,all                               ! ʷ
fini 
!*******************************************
/solu
lsel,s,loc,x,kx(8000)                   ! ѡ߽߶
lsel,a,,,37,39,2
lsel,a,,,12
SFL,all,CONV,Hair,,Tair                 ! ʩӿ߽

lsel,s,,,6                              ! ѡ϶߽߶
lsel,a,,,16,17,1
SFL,all,CONV,Hyt,,Tyt                   ! ʩ϶߽
alls
solve                                   !          
finish
!*****************************************
/POST1
PLNSOL,TEMP                             ! ʾ¶ȳֲ
asel,s,,,1                              ! ѡṹ
esla                                    ! ѡṹԪ
PLNSOL,TEMP                             ! ʾṹ¶ȳֲ
FINISH
!****************еغ***********************
presi=8.83                              ! ѹ
presa=presi*(tongt_di**2)/((tongt_di+2*tongt_t)**2-tongt_di**2)  
	                                ! ѹͲ˲ƽ
mass=270000                             ! 豸
presg=4*mass*9.8/(pi*((tongt_di+2*tongt_t)**2-tongt_di**2)) 
                                        ! 豸Ͳ˲Ӧ 
!****************½ǰ********************
/PREP7
etchg,tts                               ! ȵԪתΪṹԪ
KEYOPT,1,3,1                            ! Գѡ
MP,EX,1,2.0E5                           ! ϵĵģ
MP,NUXY,1,0.3                           ! ϵĲɱ
MP,ALPX,1,1.01E-5                       ! ϵ
!*******************************************
/SOLU
LDREAD,TEMP,,,,,,RTH                    ! ǰ¶ȳֲ
TREF,Tair                               ! ο¶
dl,3,1,uy,0                             ! ȹ׶߶ʩλԼ
dl,7,1,ux,0                             ! ͷԳʩxλԼ
sfl,6,pres,presi                        ! ڱʩѹ
sfl,16,pres,presi                       ! ڱʩѹ
sfl,17,pres,presi                       ! ڱʩѹ
sfl,1,pres,presg-presa                  ! Ͳ˲ʩƽغ
asel,s,,,1                              ! ѡṹ
nsla,,1                                 ! ѡṹڵ
esla,,1                                 ! ѡṹڵ
solve                                   ! õӦֲ          
finish
!*******************************************
/POST1
plnsol,s,int,0,1                        ! Ӧǿͼ
!****************Ӧ·1--1****************
node11=node(kx(800),ky(800),kz(800))    ! ڽڵ
node12=node(kx(11),ky(11),kz(11))
path,set1,2                             ! ·Ƽ·       
ppath,1,node11                          ! ݽڵ㶨·
ppath,2,node12                          ! ݽڵ㶨·
prsect,-1                               ! г·Ի
!****************Ӧ·2--2****************
node21=node(kx(4),ky(4),kz(4))
node22=node(kx(900),ky(900),kz(900))
path,set2,2        
ppath,1,node21
ppath,2,node22
prsect,-1
!****************Ӧ·3--3****************
node31=node(kx(10),ky(10),kz(10))
node32=node(kx(7),ky(7),kz(7))
path,set3,2        
ppath,1,node31
ppath,2,node32  
prsect,fengt_di/2
fini
